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A lattice calculation of the ttN sigma term is described using dynamical staggered fermions. Preliminary results 
give a sea term comparable in magnitude to the valence term. 



1. Theoretical Discussion 

The TriV sigma term, a-^N, is defined as that 
part of the mass of the nucleon coming from the 
expectation value of the up (u) and down (d) 
quark mass terms in the QCD Hamiltonian, 



CttN — m{N\uu + dd\N), 



(1) 



where we have taken these quarks to have equal 
current mass {— to). Other contributions to the 
nucleon mass come from the chromo-electric and 
chromo-magnetic gluon pieces and the sea terms 
due to the s quarks. Experimentally this matrix 
element has been measured from low energy tt- 
N scattering, A dehcate extrapolation to the 
chiral limit 1^ gives a result for the isospin even 
amplitude of S//^ with E = (JttN, from which the 
ttA'^ sigma term may be found. The precise value 
obtained this way has been under discussion for 
many years. For orientation we shall just quote 
a range of results from later analyses of Ut^jsi « 
56MeV, [|, down to 45MeV, [|. 

To estimate valence and sea contributions to 
cr^AT, classical current algebra analyses assume 
octet dominance and make first order perturba- 
tion theory about the SUpii) flavour symmetric 
Hamiltonian. This gives 



fir AT 



m{N\uu + dd — 2ss\N)syn 
+2Tn{N\ss\N)s ym 



def 



val 



sea 



(2) 



where we have first assumed that the nucleon 
wavefunction does not change much around the 
symmetric point. We then subtract and add 
a strange component. At the symmetric point 
the u and d quarks each have equal valence and 
sea part, while the s quark matrix element only 
has a sea component. Thus in the first term 
the sea contribution cancels, justifying the def- 
initions given in eq. (^). Using first order per- 
turbation theory for the baryon mass splittings 
a^ff may be calculated to give cr^;'^ « 25MeV 
and so crj'j^ « 31 ~ 20MeV. This in turn means 
that ms{N\ss\N) w 400 ~ 250MeV, which would 
indicate a sizeable portion of the nucleon mass 
(938MeV) comes from the strange quark contri- 
bution. 

2. Measuring a-^N 

We now turn to our lattice calculation. We 
have generated configurations using dynamical 
staggered fermions on a 16^ X 24 lattice at /3 = 
5.35, TO = 0.01 (plus some larger masses), 
Practically there are several possibilities open to 
us for the evaluation of the matrix element. The 
easiest is simply to differentiate the shift operator 
S4 giving 

OMn, 



(3) 
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(the Feynman-Hellmann theorem) . Thus we need 
to measure Mn for different masses to and numer- 
ically estimate the gradient. This leads to 



a^N ~ 11.9(8)toU=o.oi ~ 0.12(1), 



(4) 
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giving (jT^fq/Mjq « 0.15, which is to be com- 
pared with the experimental result of UT^jq /Mpf ~ 
0.06 ^ 0.05. The numerical result is much larger, 
but presumably this simply indicates that we 
have used much too large a quark mass in our 
simulation (m^^^ « 35MeV, |^). At present we, 
and everybody else, cannot avoid doing this. 

We now wish to separately estimate the valence 
and sea contributions. This is technically more 
complicated and involves the evaluation of a 3- 
point correlation function, (i? = baryon) 

C{t;T) = {B{t)mxx{r)Bwm, (5) 

(W = wall source). This may be diagrammati- 
cally sketched as the sum of two terms: 




t 



(a) is the connected (or valence) term; the nucleon 
quark lines are joined with the mass insertion, (b) 
represents the disconnected (or sea) term. 

For the connected piece we have fixed t to be 
8 or 9 and then evaluated C(t;T) as a function 
of r. The appropriate fit function is, with /x^ = 
^aexp(— A/q), (A, = ^1 denoting the parity 
partner to the nucleon, N, = +l)j 

C(i;T)« ^a/3(a|mxxl/3K~>^, (6) 

(|T ^ t ^ T ^ 0). Aa/3 and Ma are known from 
2-point correlation functions. We see that when 
a = P = N we have the matrix element that we 
require: {N\mxx\N)- However there are other 
terms which complicate the fit: (A|mxx|A) and 
the cross terms (A|toxxI-^) = MaA^7v^(-^|"^XX|^) 
which are responsible for oscillations in the result. 
To disentangle the wanted result from (A|toxx|A) 
we need to make a joint fit to the t ~ 8 and i = 9 



results. At present we have not done this, but just 
checked that separate fits give consistent results. 
In Fig. 1^ we show preliminary results for t = 8. 
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Figure 1. The ratio of the 3-point connected cor- 
relation function to the 2-point correlation func- 
tion against r, with a fit from 2 to 6. (The dotted 
lines just extend the fit to neighbouring points.) 

(Up until now we have only evaluated about a 
quarter of our available configurations.) We find 

0.08(2), (7) 

with about the same result for the A matrix ele- 
ment. The cross term is smaller, roughly 0.02. 

Finally we have attempted to estimate the 
disconnected term using a stochastic estimator, 
0. (100 sets of Gaussian random numbers were 
used.) One can improve the statistics by sum- 
ming the 3-point correlation function over r; this 
is then equivalent to a differentiation of the 2- 
point function with respect to m. This gives 

"YC{t;T)^t ^ Aaa{a\mxx\a)^li, (8) 

T a=N,A 

(ir > t » 0). In Fig. I we show the discon- 
nected 3-point correlation function. As expected 
the quality of the data is poor - optimistically we 
see a slope. A simple linear fit gives 

« 0.05(4). (9) 
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Figure 2. The ratio for the disconnected corre- 
lation function. To improve the signal we have 
averaged over the points t and T — t. 

3. Discussion 

We see that roughly eq. is consistent with 
eqs. (0^) and that a-^jq ja""^^ k, 1.5, which is to be 
compared with the experimental value of about 
2.2 ^ 1.8. Although wc can draw no firm conclu- 
sions at present our result tentatively indicates 
that the valence term is slightly larger than the 
sea term. 

Comparing our results with other work ob- 
tained using dynamical fermions, |^ uses 2 
flavours of Wilson fermions and finds for the ra- 
tio 2 ~ 3, which indicates a somewhat larger sea 
component. On the other hand, ||] using results 
from W\ for 2 staggered flavours finds 1.5 ^ 2.0, 
while [ |ll[ has ~ 2.0. These, like our result, seem 
to be lower than for the Wilson fermion case. 

We would also like to emphasise that although 
lattices can give a first principle calculation, at 
present one is not able to do this. Technically 
the fit formulae that we employ, eqs. ||) are 
true only for the complete correlation function. 
The best way to circumvent this problem and 
others is to make simulations at different strange 
quark masses; differentiation as in eq. (|) would 
then give directly ms{N\ss\N) ^ the strange con- 
tent of the nucleon. However this calculation is 



not feasible at the present time. This should be 
the ultimate goal of lattice simulations, as other 
recent theoretical results, have hinted that 
perhaps the strange quark content of the nucleon 
is not as large as supposed, previous results be- 
ing explained by a combination of factors, such 
as E 7^ fTjrJV and higher order corrections to first 
order perturbation theory. (Indeed there are al- 
ready tantalising lattice indications, [ p^lp^ ].) 

In conclusion we would just like to say that 
lattice results at present are generally in qualita- 
tive agreement with other theoretical and exper- 
imental results. However much improvement is 
required to be able to make quantitative predic- 
tions. 
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